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APPROXIMATE FORMULAE FOR PRICING ZERO-COUPON 
BONDS AND THEIR ASYMPTOTIC ANALYSIS 



Abstract. We analyze analytic approximation formulae for pricing zero- 
coupon bonds in the case when the short-term interest rate is driven by a 
one-factor mean-reverting process with a volatility nonlinearly depending on 
the interest rate itself. We derive the order of accuracy of the analytical ap- 
proximation due to Choi and Wirjanto. We furthemore give an explicit formula 
for a higher order approximation and we test both approximations numerically 
for a class of one-factor interest rate models. 

Key Words. One factor interest rate model, Cox-IngersoU-Ross model, bond 
price, analytical approximation formula, experimental order of convergence. 



1. Introduction 

Term structure models give the dependence of time to maturity of a discount 
bond and its present price. One-factor models are often formulated in terms of 
a stochastic differential equation for the instantaneous interest rate (short rate). 
In the theory of nonarbitrage term structure models the bond prices (yielding the 
interest rates) are given by a solution to a parabolic partial differential equation. 
The stochastic differential equation for the short rate is specified either under a 
real (observed) probability measure or risk-neutral one. A risk-neutral measure is 
an equivalent measure such that the derivative prices (bond prices in particular) 
can be computed as expected values. If the short rate process is considered with a 
real probability measure, a function A describing the so-called market price of risk 
has to be provided. The volatility part of the process is the same for both real and 
risk-neutral specification of the process. The changes in the drift term depend on 
the so called market price of risk function A. 

It is often assumed that the short rate evolves according to the following mean 
reverting stochastic differential equation 



where cr > 0, 7 > 0, a > 0, /3 are given parameters. In particular, it includes the 
well known Vasicek model (7 = 0) and Cox-IngersoU-Ross model (7 = 1/2) (c.f. 
Vasicek f?) and Cox, IngersoU and Ross (^). For those particular choices of 7 
closed form solutions of the bond pricing PDE ([2]) are known. Assuming a suitable 
form of the market price of risk it turns out that both the real and risk neutral 
processes for the short rate have the form ([1]). More details concerning the term 
structure modeling can be found in Kwok U). 
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dr — (a + Pr)dt + ar'^dw 
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Using US Treasury Bills data (June 1964 - December 1989), the real probability 
model ([!]) and generalized method of moments Chan et al. ^) estimated the pa- 
rameter 7 at the value 1.499. This is considered to be an important contribution, 
as it drew attention to a more realistic form of the short rate volatility (compared 
to Vasicek or CIR models). Using the same US Treasury Bills data, Nowman in (0) 
estimated 7 = 1.361 by means of Gaussian methodology. It should be noted that 
these estimations of 7 are beyond values 7 = or 7 = ^ for which the closed form 
solution of the bond prices is known in an explicit form. In (0) a model with interest 
rates from eight countries using generalized method of moments and quasi maxi- 
mum likelihood method has been estimated. They tested the restrictions imposed 
by Vasicek and CIR models using the J-statistics in the generalized method of mo- 
ments and likelihood ratio statistics in the quasi maximum likelihood method. In 
all tested cases except of one, the restrictions 7 = or 7 = ^ were rejected. Hence, 
the study of the bond prices for values of 7 different from and 1/2 can be justified 
by empirical results. However, in these cases no closed form expression for bond 
prices is known. An approximate analytical solution was suggested in ([l|) which 
could make the models with general 7 > to be more widely used. In this paper, 
we analyze the analytical approximation by Choi and Wirjanto (|1() and derive its 
accuracy order. Furthemore, by adding extra terms to it we derive an improved, 
higher order approximation of the bond prices. 

The paper is organized as follows. In the second section, we derive the order of 
approximation of the analytical approximative solution from 111). We derive a new, 
higher order accurate approximation. In the third section, we compare the two 
approximations with a known closed form solution from the CIR model (7 = ^). 
In Appendix we provide a proof of uniqueness of a solution of a partial differential 
equation for bond pricing for the parameter range 5 < 7 < |. 

2. Accuracy of the analytic approximation formula for the bond price in 
the one-factor interest rate model 

In ([l|) the authors proposed an approximate analytical formula for the bond 
price in a one-factor interest rate model. They considered a model having a form 
(m under the risk-neutral measure. It corresponds to the real measure process: 

dr = (a + /3r + X{t, r)ar'^) dt + ar'^dw 

where A(i, r) is the so called market price of risk. For a general market price of risk 
function A(t, r), the price P of a zero-coupon bond can be obtained from a solution 
to the following partial differential equation: 

(2) - drP + ]^<T^r^^dlP +{a + f]r)drP - rP = 0, r > 0, r e (0, T) 

satisfying the initial condition P(0, r) = 1 for all r > (see e.g. 13, Chapter 7)). 

Definition 1. By a complete solution to (0) we mean a function P — P{t, r) having 
continuous partial derivatives drP, drP, d^.P on Qt = [0, 00) x (0,T), satisfying 
equation ^ on Qt, the initial condition for r £ [0,oo) and fulfilling the following 
growth conditions: |P(t, r)| < Afe~""' and |Pr(T, r)\ < M for any r > 0, i G (0, T), 
where M, m,S > are constants. 

It is worth to note that comparison of approximate and exact solutions is mean- 
ingful only if the uniqueness of the exact solution is guaranteed. The next theorem 
gives us the uniqueness of a solution to ([2]) satisfying Definition [TJ In order not to 
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interrupt the discussion on approximate formulae for a solution to ([2]) a PDE based 
proof of the uniqueness of the exact solution is postponed to Appendix. 

Theorem 1. Assume 5<7<f or = ^ and 2a > . Then there exists a 
unique complete solution to 

Now let us state the main result on approximation of a solution to ^ due to 
Choi and Wirjanto (^1). They proposed the following approximation P°'P for the 
exact solution P'^^: 

Theorem 2. (1, Theorem 2) The approximate analytical solution P°'P is given by 

a 



lnP'^PiT,r) = -rB + -{t ~ B) + {r'-> + qr) — 



a2 



i?^(2/3r-l)-2i?(2r-|)+2r2-|: 



where q{r) = 7(27 - l)a^r^^^''-^^ + 2'^r^'<-^{a + (3r) and B{t) = (e'^^ - l)//3. 

Derivation of the formula ^ is based on calculating the price as an expected 
value under a risk neutral measure. The tree property of conditional expectation 
was used and the integral appearing in the exact price was approximated to obtain 
a closed form approximation. 

Authors furthermore showed that such an approximation coincides with the exact 
solution in the case of the Vasicek model. Moreover, they compared the above 
approximation with the exact solution of the CIR model which is also known in a 
closed form (c.f. (Q)). Graphical and tabular description of the relative error in the 
bond prices has been also provided in (jl|). 

The main purpose of this paper is to derive the order of accuracy of the approx- 
imation formula ([3]) by estimating the difference In P°'P — In P'^^ of logarithms of 
approximative and exact solutions of the bond valuation equation ([2]). Then, we 
give an approximation formula of higher order and we analyze its order of conver- 
gence analytically and numerically. 

2.1. Error estimates for the approximate analytical solution. In this part 
we derive the order of accuracy for the approximation derived by Choi and Wirjanto 

a. 

Theorem 3. Let P^^ be the approximative solution given by and P"^^ be the 
exact bond price given as a unique complete solution to 0). Then 

\nP''P{T,r) - In P'=^(T,r) = C5(r)r^ + o(r^) 

as T — > 0^ where 

C5(r) = _-L^^2(7-2)^2 ^2a^-l + 2^y + ip^jr* - Sr^+^V^ 

(4) +2/3(1 - 57 + 672)r2(i+->')(72 + a'^r^-r(2-/ - 1)2(47 - 3) 

+2ar (/3(-l + 4-f)r^ + {2j - 1)(37 - 2^''^^)] . 

The convergence is uniform w. r. to r on compact subintervals [ri,r2] CC (0,cxd). 

Remark 1. The function c^{r) remains bounded as r for the case of the CIR 

model in which 7 = 1/2. More precisely, limr^oC5(?') = —j^afi. If 1/2 < 7 < 1, 
then C5(r) becomes singular, C5(r) = O (r^^''^^^) as r ^ 0+. 
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Proof: Recall that the exact bond price P"^^ (r, r) for the model ([T]) is given by a 
solution of the PDE Let us define the following auxiliary function: /"^^(r, r) = 
lnP'=^(r,r) . Clearly, drP^'' = P^'-'drf^'-' ,drP^'-' = P'^^drf^ and 



q2 pex pex 

follows: 



i2 fex 



Hence the PDE for the function f^'-^ reads as 



(5) 



drr 



1 



2 2'Y 

a r ' 



i2 ^ea: 



+ (a + pr)drr'' - r = 0. 



Substitution of /"p = lnP°P into equation (l5|) yields a nontrivial right-hand side 
hir^r) for the equation for the approximative solution f""^: 



(6) - d^r + ^^^r^^ 



Q2jap 



{a + I3r)drn ~r = h{T, r). 



If we insert the approximate solution into Q then, after long but straightforward 
calculations based on expansion of all terms into a Taylor series in r we obtain: 

(7) h{T,r)^ki{r)T^ + h{r)T^ +o{t^) 

where fc4 and fcs are given by 
1 



24 



(8) 



+2/3(1 - 57 + 672)r2(i+'^)cr2 0.4^47 ^ _ 23^2 _^ -^g^a^ 
+2ar + 47)r2 + (2 - 77 + Q-f'^y'^'' a"^)] , 



6a^/3 (-1 + 27) + UP^-^r'^ - 10(1 - 27)^ri+''''a'* 



2^^2(-2+7) 

120 

+6/3^^(1-57 + 672) r2(i+T) 



(9) 



+/3r2^cr2 |^_j^o (5 + 27) + 3(1 - 27)^ (-3 + 47) r^Tci^^ 
+2Q!r ^3/3^ (-1 + 47) + 3/3 (2 - 77 + 67^) r^^cr^ 
- 5 (-1 + 27) r^+^V^' 



Let us consider a function g(T, r) = jap _ jex _ (drg)"^ = {drf°'P)'^ — {drf^^)'^ 



2drf'^^drg we have 
~drg 



^ 2 2'Y 

—crr^ 
2 



idrgf + {d^^g)\+ia + Pr)drg 
1 



- -drr 



2 27 



1 



2 27 
-CT r ' 



{drrf+{d^rr)]+io^ + f3r)drf 



-a^^-^ Org ■ 

It follows from ^ and ([5]) that the function g satisfies the following PDE: we obtain 
a PDE for the function g: 

- drg 



1 2 27 



(drgY +d'rg\+{a + pr)drg 
(10) = hiT,r)-a'r'''{drDidrg), 

where h{T,r) satisfies ([7]). Let us expand the solution of (|10p into a Taylor sc- 
ries with respect to r with coefficients depending on r. We obtain g(r, r) = 
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X^i^o "^j(^)'''* ~ '^iZuj the first nonzero term in the expansion is C(^(r)T". 

Then drg = a;c„(r)r'^-i + o(t"-1) and /i(r, r) = k4{r)T^ + o(r'*) as r ^ 0+. Here 
the term fc4(r) is given by ([5]). The remaining terms in ((T]) are of the order o{t'^~^) 
as r ^ 0^. Hence — wc„(r) = A:4 (r)T* from which we deduce, for uj — 5, c^{r) = 
-■ifc4(r). It means that g(r,r) = In P'^P(T,r) - lnP'=^(r,r) ^ -^k4{r)T^ + o{t^) 
which completes the proof. <^ 

Corollary 1. Theorem \^ enables us to compute error in yield curves which are 
given by R{T,r) = _]l}JlkL£l ^jkI relative error in bond prices. 

(1) The error in yield curves can be expressed as 

R^PiT, r) - R^^{t, r) = -C5{r)T^ + o(t^) as r ^ 0+; 

(2) The relative erro)0 o/ P is given by 

P''nr,r)-P^^{r,r) ,,5^,5, n+ 
^ = -C5(r)T + o(r ) as T ^ 0+. 

The convergence is uniform w. r. to r on compact subintervals [ri,r2\ CC (0,cxi). 

Proof: The first corollary follows from the formula for calculating yield curves. 
To prove the second statement we note that Theorem |3] gives In — In P'^^ = 
C5{r)T^ +o{t^). Hence P^p/P^^ = gC5(r)r=+o(r=) ^ ij^^^(^^y5 _^_g(^^5-^ therefore 
P^=^^-c,{r)r^+o{T^). 

Remark 2. For the CIR model with j — 1/2 the term fc4(r) defined in (0) can be 
simplified to ^cr^ [a/3 + r(/3^ — 4ct'^)] and hence 

\nP^%{T,r) - lnPa7fi(r,r) = -^a' [aP + r{l3^ - V)] + o{t^) 

as T — !• 0+ uniformly w. r. to r on compact subintervals [ri,r2\ CC [0,cx)). 

2.2. Improved higher order approximation formula. It follows from ([3]) that 
the term In P"P(r, r) — cc,(r)T^ is the higher order accurate approximation of InP*^^ 
when compared to the original approximation In P''P(T,r) from ([l|). Furthemore, 
we show, that it is even possible to compute 0{t^) term and to obtain a new 
approximation In P'^p^ (^^^ j,-) ^^q\^ thai the difference In P'^p^ (r, r) — In P*^^ (r, r) is 
o(r^) for small values of r > 0. 

Let P'^^ be the exact bond price in the model ([T]). Let us define an improved 
approximation P^p^ ^yy the formula 

(11) lnP"P2(r,r) = In P"P(r,r) - C5(r)r5 - CQ{r)T^ 

where In P°p is given by ([3]) , cs (r) is given by (U) in Theorem 1 and 

c,{r) = 1 QaV2^4'(r) + {a + Pr)c',{r) - k,{r)^ 

where c'^ and Cg stand for the first and second derivative of c^{r) w. r. to r and 
is defined in ([9]). 

Theorem 4. The difference between the higher order approximation In P^p^ given 
by ill]] and the exact solution InP*^^ satisfies In P''P^(r,r) — In P'^^(T,r) = o(r^) 
as T — !■ 0+. The convergence is uniform w. r. to r on compact subintervals 
[ri,r2] CC (0,oo). 



^This is referred to as the relative mispricing in 
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Proof: We have to prove that g{T,r) = C5(r)T^ +cg(r)r^ +o(r^) where C5 and Cg 
are given above. We aheady know the form of the coefficient C5 = c^if)- Consider 
the foUowing Taylor series expansions: 



00 

i— 5 i—A i—1 

The absolute term Iq is zero because /'^^(O, r) ~ In P^^(0, r) = In f = for all r > 0. 
Substituting power series into equation (jlOp and comparing coefficients of the order 

enables us to derive the identity: — 6c6(r)+i(T^r^'''c5 (r)+(a+/3r)c5(7')— fc5(r) = 
and hence ce{r) — ^ [^(T^r^''c'l{r) + {a + (ir)c'^{r) — k^{r)) The term k^{r) given 
by ([9]) is obtained by computing the expansion oi h. <0> 

The order of relative error of bond prices and order of error of interest rates for 
the new higher order approximation can be derived similarly as in Corollary [1] 

Remark 3. It is not obvious how to obtain the next higher order terms of expansion 
because the equations contain unknown coefficients hir), i>l, of logarithm of the 
exact solution which is not known explicitly. 

Remark 4. In the case of the CIR model we have 

and so c^^^ir) = ^ {~2af3^ + 17/3crV - 2p^r + 2acr^) . Hence 



,2 



InP^f^^lnP^^,^ + -_ (a/? + r(/?2 - 4a^)) 

2 

{-2al3^ + UPa'^r - 2/3V + 2aa^) 

The theorem yields In P^j^(r,r) — lnPpj^(T,r) = o(r^). By computing the expan- 
sions of both exact and this approximative solutions we finally obtain 

2 / 

lnP^f^(r, r) - lnPa-,;,(T, r) ~ ^ [lla(3' + 11/3^ - SAa^a' 

-180l3\a^ + 34rcr'^^T'^ + o{t'') as r ^ 0+. 

2.3. Comparison of approximations to the exact solution for the CIR 
model. In this section we present a comparison of the original and improved ap- 
proximations in the case of the CIR model where the exact solution is known. We 
use the parameter values from (1), i.e. a = 0.00315, /3 = —0.0555 and a = 0.0894. 

In Table [T] we show Loo and L2 — norms with respect to r of the difference 
\-aP^P - InP^^ and In P^p^ _ inP'=^ where we considered r e [0,0.15]. Maxi- 
mum value considered 0.15 means 15 percent interest rate, which should be suf- 
ficient for practical use. We also compute the experimental order of convergence 
(EOC) in these norms. Recall that the experimental order of convergence gives 
an approximation of the exponent a of expected power law estimate for the error 
II lnP''P(r, .) - In P'=^(t, .)|| = 0(t") as t ^ 0+. The EOC, is given by a ratio 

EOQ = Hern/ern.+i) ^^^^^ ^^^^ ^ || i^pap^^^^ ^ _ In P-(t„ .)||p . 
ln[Ti/Ti+i) 

In Table [5] and Figure [1] we show the L2 — error of the difference between the 
original and improved approximations for larger values of r. It turned out that the 
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Table 1 . The Loo and L2 — errors for the original In Pcjn and 
improved InP^j^ approximations 



r 


II lnP"P - lnP"||oo 


EOC 


II lnP°f^ - lnP===||oo 


EOC 


1 

0.75 
0.5 
0.25 


2.774 X 10" ' 
6.717 X 10"** 
9.023 X 10"^ 
2.876 X 10"^° 


4.930 
4.951 
4.972 


4.682 X lO"-'" 
6.181 X 10"" 
3.576 X 10"^^ 
2.786 X 10"" 


7.039 
7.029 
7.004 




T 


II InP"" - In II 2 


EOC 


II InP""^ - InP'^lla 


EOC 


1 

0.75 
0.5 
0.25 


6.345 X 10"" 
1.535 xlO"** 
2.061 XIO"" 
6.563 xlO"^^ 


4.933 
4.953 
4.973 


9.828 X 10-" 
1.296 X 10"" 
7.492 X 10"^'^ 
5.805 X 10"^^ 


7.042 
7.031 
7.012 



Table 2. The L2 — error with respect to r for large values of r. 



r 


1 


2 


3 


4 


5 


II lnP"P - lnP"||2 
II InP""^ -InP-'^^lla 


6.345 X 10"" 
9.828 X 10"" 


1.877 X 10"'= 
1.314 X 10"" 


1.314 X 10"^ 
2.329 X 10"^ 


5.093 X lO"'' 
1.799 X 10"'^ 


1.427 X 10"* 
8.798 X 10"*^ 




r 


6 


7 


8 


9 


10 


[| lnP"P - lnP="||2 
||lnP°p2 _lnP=^||2 


3.255 X 10"* 
3.217 X 10"^ 


6.441 X 10"* 
9.618 X 10"^ 


1.148 X 10"-* 
2.479 X 10"* 


1.890 X 10"-^ 
5.705 X 10"* 


2.921 X 10"-* 
1.200 X lO"-'* 




2 4 6 8 10 



Figure 1. The error ||lnP"P(r, .) - lnP''^(r, .)||2 for the origi- 
nal approximation (dashed line) and the new approximation (solid 
line). Horizontal axis is time to maturity r. 

higher order approximation P^p^ gjyes about twice better approximation of bond 
prices in the long time horizon up to 10 years. 

2.4. Comparison of approximate and numerical solutions. In Table [3] we 
present a comparison of the original approximation formula with a numerical so- 
lution p""™. The numerical solution was obtained using a finite volume method. 
We used 10^ spatial and 4T0'' time discretization grid points in the computational 
domain r S [0,1], r G [0,0.5] in order to achieve the L2 — errors less than 10^^^ 
between exact solution for the CIR model and the numerical solution. The differ- 
ence 0(10^^^) between the numerical and approximate solutions is therefore of the 
same order of accuracy as the numerical scheme and hence it was not reasonable 
to compute EOC in this case. 
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Table 3. Norms of the difference In P'^p(t, .) - lnP""™(T, .) for 
several values of r and 7. 





7 = 0.5 


7 = 0.75 


r 


Loo norm 


L2 norm 


Loo norm 


L2 norm 


1 

0.75 
0.5 
0.25 


2.771 X 10"' 
6.694 X 10"** 
8.854 X 10"'' 
3.400 X 10"^" 


8.967 X 10"** 
2.165 X 10"** 
2.867 X lO"'' 
7.236 X 10"" 


5.576 X 10"** 
1.691 X 10"** 
1.411 X lO"'^ 
6.963 X 10"" 


1.429 X 10"** 
3.429 X 10"" 
4.656 X 10"^° 
9.542 X 10"" 






7 = 1.00 


7 = 1.32 


r 


Loo norm 


L2 norm 


Loo norm 


L2 norm 


1 

0.75 
0.5 
0.25 


5.798 X 10"" 
1.216 X lO"'' 
9.071 X 10"^" 
6.154 X 10"^" 


1.296 X 10"" 
2.838 X 10"^" 
7.488 X 10"" 
5.663 X 10"" 


2.664 X 10~" 
1.406 X 10"^ 
1.113 X 10"^ 
7.860 X 10"^" 


5.536 X 10"'" 
2.352 X 10"^'' 
1.413 X 10"^" 
8.524 X 10"" 



3. Conclusions 

We analyzed qualitative properties of the approximation formula for pricing zero 
coupon bonds due to Choi and Wirjanto (P). We furthermore proposed a higher 
order approximation formula for pricing zero coupon bonds. We derived the order 
accuracy for both approximations and we test them numerically. The improved 
approximation is more accurate for a reasonable range of time horizons. 
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Appendix A. Uniqueness of a solution to zero coupon bond PDE 

In this section, we give a proof of Theorem 1 . Our aim is to prove the inequality 

(12) — / r'^P'^dr < K / r'^P^dr 

dr J a Jo 

to be satisfied by any solution of ^ with some constants K and u; > 0. It implies 
the uniqueness of a solution to the PDE Indeed, if Pi and P2 are two solutions of 

(12) with the same initial condition P{0, r) — 1. Then P = Pi — P2 is also a solution 
to (HI) with P(0,r) = 0. Let us define a function j/(t) = r'^ P'^{T,r)dr. Then 

the inequality (I12|) means ^^^7-^ < Ky{T) for r > 0. It implies: ^ {s^^^y{T)) = 
-Ke-^^y{T) + e-^^^^ < 0. Since y{0) = and ?/(t) > 0, it follows that ?/(r) = 
for all T. Thereof P{t, r) — for all r > 0, r > and hence Pi = P2 as claimed. 
Now let us derive inequality (|12|) . Multiplying the equation by r'^P, where 
> and 27 + cj - 1 > using the identity r'^P'^dr = r'^PdrPdr, and 

integrating with respect to r from to infinity we obtaiiQ 

(13) _ / r-'ps^^ / r'^-t+'^df.PP + / {a + f3r)r'^drPP - / r'^+^P^. 
2 dr Jo 2 Jq 7o Jo 



'In what follows, we shall omit the differential dr from the notation 
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We use the notation P' = drP-, P" = d^P. Firstly, we use integration by parts 
for the following integrals from the above equation: 



/ r^-r+'^P"P = -{2j + Lu) / r^-'+^-ipP'- / 
Jo Jo Jo 



'{P'f 



2 

where we have used the identity Jq°° r'^+ip'p = — /g°° j-'^+i-^p'^ valid for any 
w, ^ > and a function P satisfying the decay estimates from Definition [1] Substi- 
tuting this to (1131), we end up with the identity 

r-P' = ^-{2^ + u){2^ + U-l)j^ ,2,+.-2p2_^y^ ,2,+.(p,)2 

(14) - — r r^-^P^-tmi r r'^P^- r r^+^p\ 

2 Jo 2 Jo Jo 

Case 1: 7 = i and 2a> a^. We recall that the condition 2a > in the case of 
CIR model (7 = i) is very well understood as it almost surely guarantees the strict 
positivity of the stochastic processes r = rt satisfying the stochastic differential 
equation: d?- — [a + /3r) dt + a^/rdw (see e.g. 13)). 

Subcase la: 2a > a'^. We use the equality ([14|) with 7=1/2 and w = 2" _ 1 > 

to obtain the desired inequality with K — {u; + 

Subcase lb: 2a = a^. Using identity (|14p with w = (or simply by multiplying 
the PDE with P and integrating over (0, 00)) we obtain the inequality with 
K = p. 

Case 2: 76 (^j l)- We use equation with uj — 2 and estimate the integral 
J- ^27 p2 by 

using Holder's inequality: 

poo poo / poo \ 27 — 1 / poo \ 2 — 27 

j r^^P^ = J (r4^-2p47-2) 1^^2-27 p4-47)^ < (^J ^2 p2^^ j^j ^p2 j 

It follows from the Young's inequality ah < -\aP + -e^V^ for p,q > 1 such that 

1 -I- i = 1 and any e > we get 



r 



/-oo 



< (27 - 1) - / r^P^ + (2 - 27)£^ / rP 



Again using (fH)) with lo ^ 2 and the above estimate we obtain 

' r^P^ < ^(7 + l)(27 + l) / r^-'P^^a rP^ - ^ I 



At Jq 2 Jq Jq 2 

< K r^P^+(a^{j+l){2-f + l){l--f)e'^-a^J rP^ . 

where K = ^{^i + \){2-i + l)(2-i - I) {\)^ - ^. By choosing e > sufficiently 
small such that ct^(7 + 1)(27 + 1)(1 — 7)22-2", — a < 0, we finally obtain the desired 
inequality /J" r^P^ < K r^P\ 

Case 3: 7=1. We again use the equation (fT4)) with oj ~ 2. we obtain (fT2)) 
withX = 3(2cr2_/5). 

Case 4: 7 g (l, |). Similarly as in the case ^ < 7 < 1 we make use of the 
Holder's inequality integral estimation: 

x; poo / />oo \ 3 — 27 / ^oo \ 27 — 2 

j.2fp2^j (r6-47p6-47) (^67-6p47-4^ < M ^2p2 \ M ^3p2 \ 
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and, by Young's inequality, we obtain, for any e > 0, 

^27 p2 < _ 27) Q j ' f r^p^ + (27 - 2)£^ / r^P^ 



1 r\ /"CC 2 /"OO q/3 poo nOC 



By (|T4| with = 2 we have 

< Kj r^P^+[a^{"f+l){2"f + l){"f-l)e^-ljJ r^P^. 

where X = ^(7 + 1)(27 + 1)(3 - 2j) (i) ^ - ^. By choosing e > sufficiently 
small such that 0-^(7 + 1)(27 + 1)(7 ~ l)e'^ - 1 < we end up with the desired 
inequality if r^P^ < K t^P\ 
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